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In this article, we revisit the century-old question of the minimal set of observables needed to 
identify a quantum state: here, we replace the natural coincidences in their spectra by effective 
ones, induced by an imperfect measurement. We show that if the detection error is smaller than 
the mean level spacing, then two observables with Poisson spectra will suffice, no matter how large 
the system is. The primary target of our findings is the integrable (i.e. exactly solvable) quantum 
systems whose spectra do obey the Poisson statistics. We also consider the implications of our 
findings for classical pattern recognition techniques. 
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INTRODUCTION 


r(l) 




In quantum physics, the state of the system can 
be unambiguously determined by measuring several 
select integrals of motion (observables conserved in 
time evolution); such observables are said to form a 
complete set of commuting observables (CSCO) [i|. 
Mathematically, the spectra of the allowed values 
of the members of the set ..., 

are the discrete sequences of real numbers, 


{{/f \ ...}, {/f \ 




r(^) 




that obey the following property: for any pair of indices 
(m, 712), there exists at least one member of the CSCO 
(say for which the corresponding elements are non¬ 


degenerate (i.e. distinct, Operationally, 

for any h, the knowledge of the real-nnmber sequence 


r(^) ^ 


f r(l) r(2) 
was. 


r(^) 




is sufficient to infer what h 


In this paper we consider a situation where the indis- 
tinguishability occurs due to an insufficient accuracy of 
the detection. We assume that the spectra of the CSCO 
members are drawn from random processes, later cho¬ 
sen to be of the Poisson type. We further assume that 
for a given observable /, two of its possible values, 1^ 
and can not be resolved if they are separated by a 
distance less than the detection error, |/„j — /nal < A/; 
in this case, they are considered to be degenerate for all 
practical purposes. The principal goal of this paper is to 
determine the probability for two observables, and 
to form an CSCO (see examples depicted at Figs. [T] 
and [2]). 


Furthermore, we will concentrate on quantum systems 
with finite spectra of size N, i.e. we will assume that a 
given system can be in any of N available states, where 
N is finite. In real-world applications, N grows with the 
system size. 


This study is partially inspired by a related issue of the 
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FIG. 1: An example of two observables, and form¬ 
ing a complete set of commuting observables (CSCO). Even 
if a particular pair of measured values of one of the observ¬ 
ables is degenerate (points in ovals), i.e. indistinct given the 
measurement error, it will not be degenerate vis-a-vis another 
observable. If both observables are measured with respective 
errors and the state of the system, n, can be 

determined unambiguously. 


functional dependence between the quantum observables 
(i.e. the ability to predict the value of one observable 
having had measured the values of several others): while 
being a clear and extremely useful concept in classical 
physics, the practical significance of the notion of func¬ 
tional dependence in the quantum world is questionable 
at best. Already in 1929, von Neumann showed Q that 
for any quantum system with N states, one can con- 
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FIG. 2: Two observables that do not form a CSCO. If 
and are the only two observables that are measured, the 
states n = 1 and n = 3 are indistinct. 


struct a set of as many as iV — 1 conserved quantities 
each functionally independent from the other TV — 2; this 
is provided that any number of degeneracies is allowed. 
On the other hand, Sutherland argued Q that any con¬ 
served quantity is functionally dependent on any other 
conserved quantity with a non-degenerate spectrum. Be¬ 
ing able to determine the state of the system amounts to 
being able to predict the value of any other observable. 
Thus, from the functional dependence perspective we are 
looking for the probability that all conserved quantities 
be functionally dependent on two other, chosen before¬ 
hand. 


STATEMENT OF PROBLEM 


Consider two fV-element-long finite (ordered) se¬ 
quences of real numbers. 


'(1) r(l) 

1 M 2 5 • 

7(1) 

' ’ "I -^n 5 • 

jp)\ 

■ • Mat j 

'( 2 ) r( 2 ) 

1 M 2 5 • 

/( 2 ) 

' ’ "I -^n 5 • 

7(2) \ 

■ • Mat j 


Both sequences are assumed to be iV-event-long ran¬ 
domly reshuffled fragments of a Poisson process. Let 
us elaborate on what this assumption means. Consider 
monotonically increasing permutations of the above se¬ 


quences, 


/j(i) f(i) f(i) j(i)^ 

Ml ) • ■ • ) -'n ) • ■ • Mat J 


M+1 


> /(I) for any n 


( f(2) f(2) j{2) 

Ml 5 • • • Mn 5 • ■ • MiV J 

lull > for any n . 


According to the definition of the Poisson process, the 
probability of finding exactly k elements of the se¬ 
quence in an interval of a length I is 

Prob[k elements of in [/, I +1]] = 
k\ ’ 


for any initial position of the interval I. Here, is 

the expectation value of the interval between any two 
consecutive elements of the sequence: 

6IW = Mean[ii% - 1^^^] . 


Likewise, 


Prob[k elements of in [/, I +X]] = 
k\ 

Sm = Mean[P^l, - . 

An important particular property of the Poisson pro¬ 
cesses is that the intervals between two successive ele¬ 
ments of the sequence are mutually statistically indepen¬ 
dent, and they are distributed according to the exponen¬ 
tial law: 


Prob[lil, - >I]= 

PTob[i^^l, - 7^2) > I] = e-^/^ ■ 

Let us now return to the original sequences and 
that are random permutations of the and 

( — 

Ml M2 ) ■ • ■ Mri ) ■ • ■ Mat j — 

Random permutation[^7j^\ i^\ ..., i^\ ..., 7^^ j] 

/ r(2) t( 2) r(2) 7*^^)') — 

Ml M2 ) ■ • ■ Mn ) ■ • ■ ) -'at j ~ 

Random permutation[^7^^\ i^\ ..., ..., 7^^^] . 

It also follows from the properties of the Poisson pro¬ 
cesses that the 7^^^ and sequences can be approx¬ 
imately generated by dropping N uniformly distributed 
random numbers on respective intervals of lengths {N -1- 
l)<57(i) and {N + 1)WW. 
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For a given sequence, say call a pair of elements, 
Ini Slid In 2 , degenerate if the difference between them 
is less than a given number, called the detecting 

error. 

Assume that the sequence l 2 ^\ ■ ■ ■, In\ ■ • ■, 
is known, and a device (emulating a process of quantum 
measurement) produces a particular element of it, 

An observer is allowed to measure it in an attempt to de¬ 
termine what the index n was (the analogue of a quantum 
state index), using the sequence as a look-up table. For 
a perfect measurement device and in the absence of co¬ 
inciding elements in the sequence, this task can be easily 
accomplished. If however the accuracy of the measure¬ 
ment is limited, and the values separated by less than 
can not be distinguished, the value In'^ may hap¬ 
pen to belong to one of the degenerate pairs, {Ini , In})- 
If, further, there is no other degenerate pair In'^ belongs 
to, n can only be said to be equal to either ni or n 2 - If 
there are other relevant degenerate pairs, the uncertainty 
can be even greater. 

Likewise, degenerate pairs of elements of the second 
sequence is determined through 

where AI^^'> is the error of the detector that measures the 
observable. Now, both In^ and In'^ are produced, 
and both detectors are used. 

When both measurements are allowed, several ambi¬ 
guities in determining the index can be removed. How¬ 
ever, if for a given pair of indices, ni and 712 , both 
(Inl\ In}) and (/ni\ li}) constitute the respective de¬ 
generate pairs, the indices ni and n 2 can never be re¬ 
solved. The goal of this work is to determine the proba¬ 
bility that no above ambiguities exist: 

Proh\no (ni, 77 , 2 ) such that 

\l}} - l}}\ < A/(i) and \l}} - l}}\ < AJ(2) ] =? . 

Physically, we are interested in computing the probability 
that the observables I^^'^ and I^^'^ constitute a complete 
set of commuting observables (CSCO). 

A PRELIMINARY STUDY OF ISOLATED 
SEQUENCES 

Monotonically increasing sequence 

Consider a length-A fragment of a Poisson process, 
^/i, 12, ..., In, - , In'^ 

In+i > in for any n , 


with the mean spacing 51. In many respects, studying 
monotonically increasing sequences is easier than their 
random-order counterparts. This is fortunate since sev¬ 
eral important conclusions about the latter can be drawn 
using the former. 

Numerical engine 

As a model for a length-A fragment of a Poisson 
process with the mean spacing 51 we use A pseu¬ 
dorandom numbers uniformly distributed on an inter¬ 
val between 0 and (A -|- 1)51. When ordered in the 
order it has been generated, the sequence serves as 
a model for a random permutation of a Poisson se¬ 
quence, (/i, I 2 , ..., In, . •., In), the primary object of 
interest. When rearranged to a monotonically increas¬ 
ing sequence, a model for a Poisson sequence per se, 
/i, J 2 , , In^ , is generated. 

To test the numerical method, we compare, in Fig. |3l 
numerical and analytical predictions for the probability p 
of two subsequent elements of the monotonically increas¬ 
ing sequence to form a degenerate pair, with respect to 
a detection error AI. The analytic prediction, 

p = Prob[lor a given n, J„+i — in < AI] = 

^ _ g-A//37 ^ (2) 

51 ’ 

follows directly from one of the properties of the Poisson 
sequences O- The numerical and analytical predictions 
agree very well. 

From ([2]) it also follows that the average number of 
degenerate pairs is 

_ /V_VriO 

M=pNp « pN , (3) 

where 

Ap = A - I (4) 

is the number of pairs of indices with two neighboring 
elements. 

Clusters of degenerate pairs 

In what follows, the analytical predictions will be 
greatly simplified thanks to our ability to assume that 
degenerate pairs do not share elements, i.e. that a given 
element can either be isolated or belong to only one de¬ 
generate pair. It is clear that the probability of forming 
a degenerate cluster must be much smaller than the one 
for an isolated pair, because of the small probability of 
forming the latter. Yet, a quantitative assessment is due. 
Assuming that that the individual spacings between the 
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Detector window, AI [d/] (A^=100) 


FIG. 3: Probability for a given pair of consecutive elements 
of a Poisson-distributed spectrum to be unresolvable by a de¬ 
tector with an error AI (i.e. to be degenerate). 51 is the mean 
spacing between consecutive elements. Red dots: numerical 
model, uncorrelated, uniformly-distributed random numbers 
on an interval, subsequently rearranged in a monotonically 
increasing succession. Green line: theoretical prediction ©• 
Numerically, the spectrum consists of 100 elements; the the¬ 
oretical prediction is universal. Numerical points result from 
averaging over 1000 Monte Garlo realizations. 

subsequent values of / are statistically independent, we 
get: 

Pc = 

Prob[ioi a given n, {In — In-i > AI) and 
{In+i - In < A/) and 

(J„+2 -/„+!< A/)] = (1-p)/ . 

(5) 

Note that the above probability is the probability for the 
index n to be at the left end of a cluster of degenerate 
pairs containing more than one pair. Fig. 0] compares the 
above prediction with ab initio numerics: the agreement 
is remarkable. The overall conclusion is that for a small 
probability p of forming a degenerate pair, one can safely 
assume that all the degenerate pairs are isolated. This 
result will be heavily used in a subsequent treatment of 
the more physical random-order sequences. 

It is instructive to also compute the probability of hav¬ 
ing no degenerate clusters at all, for two reasons, (a) this 
is the first measure in this study that characterizes the 
/i, 72, ..., in, ..., In^ as a set, with no reference to 
the order of the elements. The probability of no clus¬ 
ters we will obtain will therefore be identical to the one 
for the (Ii, I 2 , ..., 7n, ■ • ■, In) sequence, randomly or¬ 
dered. On the other hand, the monotonically increasing 
sequences are conceptually simpler, and they allow for 
more intuition in designing analytical predictions, (b) 
Using this study we will show that while the spacings 
between neighboring members of the monotonic sequence 



FIG. 4: Probability for a given pair of consecutive elements 
of a Poisson-distributed spectrum to be degenerate and to be 
immediately followed on the right by more degenerate pairs, 
as a function of the detector error AI. Red dots: numerics. 
Green line: theoretical prediction It does not depend on 
the sequence length N. The rest of the data is the same as in 

Fig. El 

are statistically independent, the degenerate pairs they 
produce are not statistically independent from the point 
of view of the randomly permuted sequence. 

Even though the spacings between the consecutive el¬ 
ements of the I sequence are statistically independent, 
the appearance of a cluster with the leftmost edge at an 
index n will be preventing an appearance of another clus¬ 
ter at n -|- 1. Thus, appearances of clusters at different 
points are not statistically independent. Nevertheless, 
for a very low cluster probability, a typical distance be¬ 
tween the clusters becomes much larger than their typical 
length, and the above correlations can be neglected. Our 
estimate for the probability of having no clusters at all 
then reads: 

Tno clusters = 

Prob[ioT all n, (in+i — In > A/) or 
(^{in+i - In < A/) and 
{In - /n_i > A/) and 

{In+2 — In+1 > A/)^] ~ 

( 6 ) 

where Np is given by Q. Above, we neglected limita¬ 
tions on the cluster length closer to the right edge of the 
spectrum, since we are mainly interested in the N ^ 00 
limit. 

As we have mentioned above, the result ([6|) applies for 
both / and / sequences, being a characteristics of these 
sequences as a set, not as sequences with a particular 
order. It is tempting to reinterpret it as the probability 
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FIG. 5: The probability of having no clusters of consecutive 
degenerate pairs of elements at all, as a function of the de¬ 
tector error Al. Red dots: numerics. Green line: theoretical 
prediction ©• Blue line: a naive combinatorial hypothesis 
®. At small detector window widths, the formula involves 
factorials of large numbers and becomes numerically unreli¬ 
able. Purple line: the same as the blue line but using an 
asymptotic expression. The rest of the data is the same as in 
Fig. El 


that M pairs chosen uniformly at random from 

N(N -1) 


Mp = 


(7) 


pairs of elements of an TV-member-strong set have no 
elements in common. The number of ways to choose M 
pairs with no elements in common is 

2^N\ 

M\(N -2M)\ ■ 


set. The physical realization of our model will involve 
two conserved quantities in an integrable system. Those 
are indeed known to realize Poisson processes [1, Q , but 
only after a permutation to a monotonically increasing 
sequence. For two observables, these permutations are 
completely unrelated. Thus a physically more relevant 
model would involve involve sequences that are repre¬ 
sented by a random permutation of a Poisson sequence, 

(/i, /2, , In) = 

Random permutation[^Ji, J 2 , ..., ..., In'^] ■ 

Numerically, those are modeled by sets of real num¬ 
bers randomly independently distributed over an inter¬ 
val, with no subsequent reordering. 

To make a theoretical prediction for the probability for 
a given pair of two subsequent elements of the sequence 
to be degenerate, we first assume that the probability 
for degenerate pairs to share elements {i.e. form clusters 
of degenerate pairs) is low (see Fig. |4] and the expres¬ 
sion ([S])). In this case, one can assume that each de¬ 
generate pair constitutes a pair of consecutive elements 
in the monotonically increasing counterpart of the se¬ 
quence. The probability of the later is 2/Np (with Np 
given by ([4])). This probability must then be multiplied 
by the probability of a neighboring pair in the monotonic 
sequence to be degenerate, p. We get 


Prandomized — 

Prob[ior a given n, In+l In ^ ^7] — 


V—>00 2 


( 10 ) 


N, 


-P 


N^- 


The number of ways to chose any set of M pairs is 

Mp\ 

M\{Mp-M)\ ■ 

Dividing the former by the latter we get the following 
expression; 

Fno clusters, conjecture — 

p—^0, Al—>00,A^—>const 9 2 AT 

Ri e ‘‘P , (9) 

where M is given by (©I. Comparing (|5]) to (O one 
notices that the former overestimates the probability of 
cluster formation by a factor of two. An interpretation of 
this simple relationship requires future research. Figure[5] 
illustrates this difference. 


2^N\{Np-M)\ 
Np\{N -2M)\ 


Fig. ini demonstrates that this prediction is indeed correct. 

We are now in the position to assess the ability of a 
single observable I to serve as a CSCO. The probability 
of that is given by the probability of having no degener¬ 
ate pairs at all. Since this probability is a property of the 
sequence as a set, it can be estimated using the monoton¬ 
ically increasing counterpart. There, the probability of 
not having degenerate pairs is the probability that none 
of the — 1 pairs of consecutive indices are degenerate. 
We get 


tPcsco, 1 = 

Pro6[for no n, In+l In ^ ^I] — 

(1 p->0, Al-s-oo^pAl^COnst ^-pAT 


( 11 ) 


Implications for a randomly permuted counterpart 

The principal goal of our paper is to assess the ability 
of two quantum conserved quantities to serve as a CSCO 


Fig. [7] shows that this probability approaches unity only 
for a detector error as low as the inverse sequence length 
N. Furthermore, the expression (HU shows that even if 
for a given N the given observable does form a CSCO, 
for larger N this property dissappears. 
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Detector window. A/ [5/] (A^=100) 

FIG. 6: Probability for a given pair of consecutive elements of 
a Poisson distributed spectrum, subsequently randomly per¬ 
muted, to be unresolvable by a detector with an error A7 
(i.e. to be degenerate). 51 is the mean spacing between con¬ 
secutive elements before the random permutation. Red dots: 
numerical model, uncorrelated uniformly distributed random 
numbers on an interval, with no further rearrangement. Green 
line: theoretical prediction CQl). Numerically, the spectrum 
consists of 100 elements. Numerical points result from an 
averaging over 1000 Monte Garlo realizations. 
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0.01 0.02 0.05 0.10 0.20 0.50 1.00 2.00 
Detector window. A/ [57] (A^=100) 

FIG. 7: Probability of having no degenerate pairs of ele¬ 
ments at all, in a Poisson distributed spectrum (regardless of 
whether reshuffled or not), as a function of the detector error 
A7. Physically, this probability corresponds to the probabil¬ 
ity for a given observable to form a complete set of commuting 
observables (CSCO). Red dots: numerics. Green line: theo¬ 
retical prediction (HU. The rest of the data is the same as in 

Fig. El 



TWO SEQUENCES 


We are finally ready to address the principal question 
posed: what is the probability that two observables with 
Poisson spectra, subsequently independently randomly 
permuted, and 7*-^^ constitute a complete set of com¬ 
muting observables (CSCO). Mathematically, the corre¬ 


sponding probability is: 

Cbsco, 2 = 

Proh[ioT no n, (7,^^^ — 7,^^^ < A7^^^) and 

(7^^i-7^)<A7(2))]= (12) 

PrandomizedPrandomized) 


N- 




where ^randomized = ^he respective probabili¬ 

ties for a given pair to be degenerate (see (fTOll ). is the 
analogue for a monotonic sequence (see ([2])), and A7(“) 
is the respective detection error. We assume the same 
spacing, <57, for both sequences. 

A naive combinatorial interpretation of this probabil¬ 
ity is the ratio between the number of ways in which 
= p(“)_/Vp pairs can be chosen from two A-element- 
long sets in such a way that no two pairs coincide and its 
unrestricted analogue: 

^bsco, 2, conjecture — 

(Afp-M(i))! (Wp-M(2))! 
(Afp)!(AAp-M(i) -M(2))! (1^) 

_2p(i)p(2) 

~ e , 

where Afp is the number of pairs of N objects and 
Np = N—1 is the number of pairs of neighboring elements 
in a sequence of length N (see ([4])). In this particular case 
the naive model generates the correct prediction. Note 
that in this case, no correlations between the degenerate 
pairs within a particular sequence are involved, thus a 
better result. 

Figure [S] allows one to compare the various predictions 
made above, for = 2A7’^^b Notice that (a) already 
at A7(i) = .2, CSCO is reached; (b) from the expres¬ 
sion m one can see that for large N, the probability 
m does not depend on the length of the spectrum N, 
and thus CSCO persist for large systems. A particular 
illustration of this phenomenon is presented in Fig. 


SUMMARY AND OUTLOOK 

In this work we show that with a large probability, 
differing from unity only by the expression in (jl3l) (see 
also Figs. [8]l9l) . two integrals of motion with Poisson dis¬ 
tributed spectra form a complete set of commuting ob¬ 
servables (CSCO); i.e. they can be used to unambigu¬ 
ously identify the state of the system, the detection error 
notwithstanding. For a given detector error, this proba¬ 
bility converges to a fixed value as the number of states in 
the spectrum increases. The above result is contrasted 
to an analogous result for a single observable. There, 
(a) only for a very small detector error can this observ¬ 
able be used as a CSCO, and (b) for a fixed detector 
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FIG. 8: Probability for two Poisson-distributed spectra, 

subsequently independently randomly reshuffled, not to have 
any degenerate pairs of elements in common. Physically, this 
probability corresponds to the probability for two given ob¬ 
servables with Poisson spectra to form a complete set of com¬ 
muting observables (CSCO). Red dots: numerics. Multicol¬ 
ored line: comprises the theoretical prediction m , a naive 
combinatorial hypothesis (O, and its large N asymptotic 
behavior, all three mutually indistinct. Note that the former 
and the latter are analytically shown to coincide. 



FIG. 9: Numerically computed probability for two Poisson- 
distributed spectra, subsequently independently randomly 
reshuffled, not to have any degenerate pairs of elements in 
common, as a function of the total length of the spectrum 
N. As in Fig. [8] the multicolored line corresponds at the 
same time to the theoretical prediction m and to the large- 
N asymptotic behavior of the combinatorial expression (113, 
which are shown to coincide. The corresponding detector win¬ 
dow sizes are =0.1 and = 0.05. 


error, the probability of a CSCO falls to zero exponen¬ 
tially as the number of states in the spectrum increases. 
Poisson spectra constitute a popular model for energy 
spectra of integrable (i.e. exactly solvable) quantum sys¬ 
tems with no true degeneracies ii- They can also be 
used to model spectra of other integrals of motion of inte¬ 
grable systems, provided they are substantially functions 
of at least two quantum numbers. Note that in generic 


non-inteCTable quantum systems, energy levels repel each 
other [1,1^. There, even a poor energy detector, with an 
error as large as the energy level spacing, would be able to 
identify a state: energy can thus serve as a CSCO. There¬ 
fore, the physical implication of our principal result can 
be formulated as follows: given a reasonably small (i.e. 
smaller than the mean spacing between levels) detector 
error, two generic integrals of motion of an integrable 
quantum system can be used as a complete set of com¬ 
muting observables, no matter how large the system is. 
On the other hand, if only one integral of motion is used, 
the appearance of unidentifiable states is unavoidable for 
large systems. 

From the functional dependence perspective, we 
showed that in poorly measured integrable systems, all 
conserved quantities are (with a close-to-unity probabil¬ 
ity) functionally dependent on any two a priori chosen 
generic conserved quantities. 

Results obtained in our article may find applications 
beyond the quantum state identification. They should 
be generally applicable in a standard pattern recognition 
setting where an unknown object must be identified indi¬ 
rectly by one, two, or more attributes. For example, us¬ 
ing our results, we can show that chances of an ambiguity 
in identifying a person who was born on a particular date 
H in a particular town Ti within a specific time interval It 
and who is currently living on a particular street S in an¬ 
other town T 2 depends neither on the length of the time 
interval It, nor on the number of the streets in T 2 , but 
solely on the number of birth per day in Ti, probability 
of a further migration to T 2 , percentage of Ti-born citi¬ 
zens in T 2 , and the average street population in T 2 ; in the 
other words, this probability depends only on the relative 
measures, the absolute measures being irrelevant. The 
value for this probability is trivially obtainable from the 
central result of this paper. More generally, for a large 
enough set of N patterns with two Poisson distributed nu¬ 
merical identifiers, each measured with a finite error, the 
probability Pcsco, 2 of having no unidentifiable patterns 
depends only on the probability p^^'i of a given value of an 
identifier to be indistinct from its neighbor above and the 
analogous probability p^‘^1 for another identifier and not 
on size of the set N. At the same time, if only a date of 
birth is known, the procedure will start frequently pro¬ 
duce ambiguous results for long enough interval of time 
considered: here, the probability of finding two people 
born in a given town on the same date within this inter¬ 
val of time approaches unity. Generally, a single Poisson 
identifier becomes unreliable for larger sets of patterns. 
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